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Abstract- In this paper we propose a learning based 
iterative numerical method for computing a root^ of a 

non-linear equation of the form f(x) = 0 in the interval 
[a,b]. Lag based method through Hermite 
interpolation modeled root discovery approach is 
developed and demonstrated in this paper. The new 
method has been tested for a series of functions 
considered by several researchers. The numerical 
experiments show that the new method is effective 
and converges faster to a root with lesser number of 
iterations and lesser number of function evaluations. 
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1. Introduction 



Intelligent searching strategies will enhance the 
performance and efficiency of an iterative numerical 
method for solving a nonlinear equation f(x)=0. 

If f is a continuous function on the interval [a,b], the 
Bisection method converges to a root of f, which 
repeatedly bisects an interval, then selects a 
subinterval in which a root must lie for further 
processing. This is called a bracket of a root and this 
method is not constrained by any other parameter 
other than sign. So, this method is relatively slow. 



Secant method uses a succession 
of roots of secant lines to better approximate a root of 
a function f(x). The Secant method does not require 
that the root remain bracketed like the Bisection 
method, and hence it does not always converge but 
reduces computational time. 

The classical Regula falsi method finds a simple 
root of the nonlinear equation f(x)=0 by repeated 
linear interpolation between the two current 
bracketing estimates. Regula falsi method combines 
features from the Bisection method and the Secant 
method, here sign as well as magnitude is used. If the 
initial end points a and b are chosen such that f(a) 
and f(b) are of opposite signs, then this assumption 
guarantees the existence of a zero of f(x) in the 
interval [a,b]. This means that the Regula falsi method 
always converges. In Regula falsi method, at every 
iteration one point is dropped, the information 
available at points which are dropped recently can be 
used for approximating the function with better 
accuracy. This enhances the knowledge in the current 
support about the function. 

Alefeld and Potra [1,2] proposed three efficient 
methods for enclosing a simple zero £ of a 

continuous function f(x) in the interval [a,b], provided 
that f(a) f(b)< 0. Starting with the initial enclosing 
interval [ai.bi] = [a,b], this method fits a quadratic 
polynomial to produce a sequence of intervals such 

that {(b n -a n )}~ =1 

£ e [a n+1 ,b„ +1 ]c[a„,bjc....c[a 1 ,fi]c[a, l b] _ 
11m (b n -a„) =0 

n— »■» 
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In a subsequent study by Frontini [6], a numerically 
comparable method through fitting a Hermite 
interpolating polynomial was suggested. Here in this 
method, by taking x 0 = a and x 1 =b as the starting 
points, with f(x 0 )f(xi)< 0, then x, +1 is taken as the root 
of Hermite polynomial P2[f,x M ,x,](x) which belongs 
to (x/_i,x,). Use of the dichotomy condition f(x^)f(x,) < 
0 ensures the existence of at least one root of P 2 [f, 
x,-i, xj(x) belonging to (x M ,x,) because of the 
continuity of P 2 [f, x M , x,](x). The author further 
suggested in a similar manner fitting a cubic Hermite 
polynomial P 3 [f, x M , xj(x) and a generalized method 
satisfying the interpolation conditions: 

4» 1 [f,a,b](a) = f (k) {a) 

P 2 « [f,a,b] (b) = f ik) (b) k=\,2...n-\^[\) 

By considering x n+1 as the root ^of Hermite 

interpolating polynomial P 2 n-i[/ r , *n-i. x n ](x) belonging 
to (x M ,x ( ), the order of convergence of the iterative 
method obtained in [6] is given as 



p = . Frontini [6] proposed that his 

2 

method works effectively only for fitting a cubic 
Hermite polynomial. 



Based on the earlier methods which emphasize as 
memory based learning [1,2] and with reinforced 
learning by making use of derivative through Hermite 
interpolation as lag based learning [6], in this paper 
we attempt to propose an iterative method as follows. 



distinct points, their function values and the values of 
their corresponding first derivatives. 

Hermite Interpolating Polynomial: The Hermite 
interpolating polynomial interpolates not only the 
function f(x), but also first order derivatives at a given 
set of data points ( Xj , f( Xj ), f'( Xj )) , i=0,1,2,...n. 

The Hermite interpolating polynomial is given by 

=^-nx-x i )i i ( Xi )][l i ( X )^ f( Xl ) 

+ i(x-xi)[iMf f ( Xi ) -K2) 

2=0 

Where 

l ^ = (x-x 0 )(x-x l )...(x-x l _ 1 )(x-x M )...(x-x n ) 
(*, -x a )(x- Xl )... (x t - x IA ) [x t - xj. . . (x, - x n ) 

We now consider fifth degree Hermite interpolating 
polynomial passing through the 

points {Xj, f(x.), f'ix^) , i=0,1,2 by putting n=2 in (2). 



Subroutine 

f _ hermite(xj , f{x,), f'(x-)) 
Xj = [a, c, b] 

f(x.)=[f(a),f(c), f(b)] 
f'(x j )=[f'(a),f'(c), f'(b)] 

x r =root of H 5 (x) 



2. The new method 



Hermite interpolation is an extension of basic 
polynomial interpolation that not only matches 
discrete information at a set of points, but also 
matches the slope (or rate of change) at those points. 

We consider x 0 = a and x-i =b as the starting points, 
with f(x 0 )f fa) < 0, a point c is taken as a<c<b. A fifth 
degree Hermite polynomial is fit for these three 



Algorithm 

1. Set a=lower bound for search, jb=upper bound for 
search such that f(a)f (b) < 0,find a point c such 
thata < c< b . Obtain corresponding f'(a), f'(b), f'(c) . 

2. Call f _hermite(x j , f{xX f'{x } )) which returns five 
roots of fifth degree Hermite polynomial. 
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3. From the above five roots select the root (say x r ) 
which lies in the interval [a,b] and for which function 
value f(x r ) is minimum. 

4. Obtain f(x r ) if f(x r ) <e , then display root= x r 

5. Compute f\x r ) 

6. If f(a) f(x r ) < 0 , then 

c = b, f{c) = f{b), f(c) = f(b) 

a = a,b = x r , f(b) = f(x r ), f'{b) = f\x r ) 

If f(b) f(x r ) < 0 , then 

c = a, f(c) = /(a), f\c) = /'(a) 

a = x r ,b = b, /(a) = f{x r ), f'(a) = f\x r ) 

Repeat step 2 until | f( x ) |<s . 

In each iteration, the information about the point to be 
dropped is stored and used as lag, (always third point 
in new method is lag) to fit fifth degree Hermite 
interpolating polynomial, hence we need only two 
function evaluations at each iteration which improves 
the efficiency of the method. 

3. The order of convergence 

Let the function f(x) defined by (n+1) points 
(x,/(x,)), i=0,1 ,2,...n, be continuous and differentiable 
(n+1) times and let f(x) be approximated by a 

polynomial <p n (x) of degree not exceeding n such 

that 

0 n (xj) = f{xj) , i=0,1,2...n, then error is given by 

E(x) = f(x)-0 n (x) 

= x n+1 (x) f[x 0 ^,..xj = M n+X (x)L 
where x n+x (x) = (x-x 0 ) (x- x x )...(x- xj 



Then by constructing a function F(x) such that 
F(x) = f(x) -<p n {x) - L7r n+1 (x) and the constant L can 
be determined so that 

L _ / (a+1) (£) 
(n + l)\ 
Hence we get 
E{x) = f{x)-<p n {x) 

= -^—* n+ l(x)f in+1) ({),x 0 <{<x n 
{n + 1) ! 

Now, we consider Hermite interpolating polynomial 

H2n+i(x), then error E(x) = f{x) - H 2n+1 (x) can be 

obtained. Here we notice that both £(x) and 

[j B+1 (i)f vanish together with the first derivatives at 

each of (n+1) points. We then form a linear 
combination of these functions 

F(x) = / (x) - H 2n+l (x) - L[;r n+l (x)f 
(2n+2) ! 

Where f(x) is assumed to have continuous 
derivatives of order (2n+2) and E, = i;{x) is in the 
interval determined by the points x,x 0 ,...,x n . 

Theorem: Let / e C 6 , at least in a neighborhood of a, 
with f (a) = 0. Let {x n } be the sequence defined as 
x n+1 = g , where £ is the root of fifth degree Hermite 
polynomial enhanced by the dichotomy procedure so 
that it converges to a with n. Let e n =x n -abe the 

error of the nth iterate. Then, e„ +1 |~ M|eJ 2 " 9196 where 
M is a positive constant. 

Proof Without loss of generality we assume 

that^_ 2 <x B _ 1 <x n . 

There exists rjs]x n _ 2 x n [ such that 
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E(a) = f(a)-H 5 [ f, x^ , x n ] (a) 
(a- x ? f (a- x , ) 2 (a- x n f i\ 

= ^ — r 6 ' 77) -> (4) 

6! 

If x n+1 is zero of ^[f.x^.x^.xjix^) in 
K_ 2 ,*J, then 

^\f^ n _ 2 x n _ l x n ]{x n+l )=Q= f{a) 

By replacing f{a) with H^[f,x n _ 2 , x n _ v x n ]{x n+l ) in 
(4), we get 

H s [f, x„_ 2 , x n _, , x n ] (x n+l )-H 5 [f,x n _ 2 , x^ ,x n ](a) 



(a - x n _ 2 ) 2 (a - x^) 2 (a - x n ) 2 
6! 



By Lagrange's theorem, there exists ge]oc,x n+1 [ such 
that 



By 



comparing (4) and (5) we get 



(« - ) = ^ (« - *„- 2 ) 2 (« - x n-i ) 2 (« - x „) 2 



Where ^ : 



/ (6) 07) 



<fe]a,^ +1 [, 



Assuming e n+1 is asymptotic to ce^,p>l, 
Expressing this in terms of E n _j , we get 



(a - * n+l )=K(a-x n _ 2 f(a- x nA f(a-x n f 

„ 2 2 2 



L fc n-l A fc n-l fc /3-l fc n-l 



This asymptotic equation is satisfied if p is the positive root of 

2 2 

=> p = — + 2 + 2p 
P 

3 2 

=> p — 2p -2p-2 = 0 

=> p = 2.9196 

then, s n+1 |~ M |e n | 2 9196 for a suitable constant M. 

It is proved in the above theorem that the new method 
converges to a root of f(x) with rate of convergence 
2.9196. Given three distinct points we obtain their 
function values and the values of first derivatives 
corresponding to these three points before starting 
the first iteration. Consecutively we need only two 
function evaluations at each iteration as lag based 
intelligent searching strategy is used which improves 
the Efficiency Index, given by 

EI= (2.9196) 7 =1.708684 -This El is better than 
secant method, Newton method, method suggested 
by Costabile [3,4] and method proposed by Frontini 
[6]. 

Table I : Comparison of El of various methods is 
shown in Table I. 



Method 


n 


P 


Vn 

El=p 


Bisection 


1 


1 


1 


Secant 


1 


1.618 


1.618 


Regula Falsi 


1 


1 


1 


Newton 


2 


2 


1.41 


Frontini 


2 


1+ ^3 


1 .6529 


New method 


2 


2.9196 


1.7087 



4. NUMERICAL EXPERIMENTS 

The results of numerical experiments are shown 
in the tables II, III and IV. 



TABLES II and III : 
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The computed results of Examples 1-4 by new method with tolerance= 1><10" 15 are given in Table II. Here Nl stands for 
number of iterations and NFE stands for number of function evaluations. Table III- displays results obtained by 
Exponential regulafalsi, Regulafalsi, Steffensen's and Newton methods(see[7]) 



s.no 


Function f (x) 


Interval for 
Search 


Root £ 


|f( ^ )| 


Nl 


NFE 


1 


Inx 


[0.5,5] 


1.00000000000000 


2.4424907e-016 


4 


13 


2 


x+ 1 - e sinx 


[1,4] 


1.69681238680976 


2.6541 269e-01 6 


3 


11 


3 


11x 11 -1 


[0.1,1] 


0.80413309750367 


6.5052130e-016 


4 


13 


4 


xe x -0.1 


[0,1] 


0.11183255915896 


4.9786564e-016 


2 


9 





Algorithm 2 (EXRF) 


Algorithm 1 (regula falsi) 


Steffensen 


Newton 




Nl 




[f(£)l 


Nl 


# If(^) 


ll 


Nl % |f(£)| 


Nl £ 


|f(£)l 


1 


7 


1.00000e+00 


0.00000e+00 


27 


1.00000e+00 


8.88178e-16 


Failure 


Divergent 




2 


11 


1.6968.e+00 


4.44089e-16 


32 


1.6968e+00 


4.44089e-16 


Failure 


Not 

convergent 




3 


8 


8.04133e-01 


4.44089e-16 


101 


8.041 33e-01 


1.25422e-13 


Divergent 


7 8.041 33e- 
01 


4.44089e-16 


4 


6 


1.1 1833e-01 


0.00000e+00 


15 


1.11833e-01 


7.4040le-16 


Failure 


Failure 
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TABLE III (see [7]) * 

'calculations given by Jinhai Chen, and Weiguo Li [7] 

TABLE IV : In this Table IV we present results through new method for some examples considered by several researchers with 
tolerance= 1x10" 15 







Interval for 










s.no 


Function f (x) 


Search 


Root 


lf( t il 


Nl 


NFE 


1 


v 2 /y 2 /^-(- Jo cinv\ Jl /1ft 
X \X /O^ v L fall IX j- \ 0 1 1 O 


[0.1,1] 


0.39942229171099 


3.1918912e-016 


3 


11 








0.13825715505682 


4.4408921 e-01 6 


3 


11 




zxe + 1 - Ze 












2 




[0,1] 


0.06931408868702 


6.9388939e-017 


4 


13 




n-o,i u,zU 


















0.03465735902085 


3.3306691e-016 


4 


13 








0.34595481584824 


4.3021 142e-01 6 


1 


7 




..2 /a „\n 

X - (1 -x) 












3 




[0,1] 


r\ O A IT A 11 TTT7 C 11 -1 

0.2451 2233375331 


9.6450b25e-01b 


4 


13 




n-o,l U,ZU 


















0.1 6492095727644 


9.3675068e-017 


6 


17 








0.51615351875793 


3.4000580e-016 


3 


11 




e (x-1)+x 












4 




[0,1] 


0.53952222690850 


5.5250943e-016 


4 


13 




n-o,1 U,ZU 


















0.55270466667922 


1.5252099e-016 


12 


29 








0.40999201798914 


2.7755576e-017 


2 


9 




x 2 +sin(x/n)-1/4 












5 




[0,1] 


0.45250914557763 


5.5511151e-017 


2 


9 




n=5,10,20 


















0.47562684859606 


4.9960036e-016 


1 


7 








0.20000000000000 


4.7184479e-016 


9 


23 




(nx-1)/((n-1)x) 












6 




[0.01,1] 


0.10000000000000 


2.498001 8e-01 6 


8 


21 




n=5,10,20 


















0.05000000000000 


2.1035805e-016 


7 


19 


7 


x 3 -e- x 


[0,1] 


0.77288295914921 


4.9960036e-016 


2 


9 


8 


xe x -1 


[0,1] 


0.56714329040978 


6.6613381e-016 


5 


15 


9 


10 x +x-4 


[0,1] 


0.53917912205280 


0.0000000e+000 


3 


11 
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5. CONCLUSION 

As learning based intelligent searching technique, we 
propose a method of fitting a fifth degree Hermite 
polynomial (for three distinct points, their respective 
function values and values of their corresponding first 
derivatives). The present method, using lag based 
searching technique with tolerance= 1*1 CF 15 has a rate 
of convergence 2.9196 and better Efficiency Index of 
1.708684. The El obtained by the present method is 
better than the El of Secant method, Newton method, 
as well as the method suggested by Costabile [3,4] and 
the method proposed by Frontini [6]. Numerical 
experiments presented in the paper show that this 
method gives better performance and it is fast 
approaching to a root of given non linear equation f(x) = 
0 with lesser number of iterations and lesser function 
evaluations and comparable to well known methods. 



[4] F. Costabile, M. I. Gualtieri and S.Serra, "An iterative 
method for the computation of the solutions of nonlinear 
equation", Calcolo vol.36, pp 17-34, 1999. 

[5] W. Gautschi, Numerical analysis. An introduction. Boston: 
Birkhauser Boston 1997. 

[6] M. Frontini, "Hermite interpolation and a new iterative 
method for the computation of the roots of non-linear 
equations", Calcolo vol.40, pp 109- 119, 2003. 

[7] Jinhai Chen, and Weiguo Li, " An improved exponential 
regula falsi methods with quadratic convergence of both 
diameter and point for solving nonlinear equations", Applied 
Numerical Mathematics vol. 57 pp 80-88, 2007. 



ACKNOWLEDGEMENT 



The authors thank ACRHEM, UoH for financial 
support. First author thanks Professor. S.P.Tewari, 
Director, ACRHEM, UoH and Professor. M. 

Sitaramaiah, emeritus Professor, MCIS, UoH, for the 
encouragement and valuable suggestions. 



REFERENCES 

[1] G. E. Alefeld and F. A. Potra, "Some efficient methods for 
enclosing simple zeros of nonlinear equations", BIT vol. 32, 
pp 334-344, 1992. 

[2] Alefeld G, F .A. Potra, Shi Y.X, "On enclosing simple 
roots of nonlinear equation", Math. Comput. Vol. 61, pp 733- 
744 (1993) 

[3] F. Costabile, M. I. Gualtieri and R. Luceri, "New iterative 
method for the computation of the solutions of nonlinear 
equations", Numer. Algorithms vol. 28, pp 87-100, 2001. 



46 



